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I.   INTRODUCTION 

In  recent  years  the  great  expansion  in  the  development 
of  semiconducting  materials  with  acceptable  thermoelectric 
characteristics  has  made  them  attractive  for,  in  special 
applications,  direct  energy  conversion  and  heat  pumping. 

Thermoelectric  power  sources  and  heat  pumping  devices 
can  offer  significant  advantages  with  respect  to  long  life, 
years  of  unattended  operation,  freedom  from  environmental 
restraints,  and  unusually  favorable  specific  power  ratios. 

It  is  true  that  even  now  thermoelectric  cooling  and 
power  sources  are  not  serious  competitors  of  powerful  heat 
pumps  and  conventional  power  apparatus ,  but  they  are 
irreplaceable  in  those  cases  where  the  overriding  require- 
ments are  small  size,  simplicity  of  construction,  quiet 
operation,  weight,  etc. 

Various  microref rigerators ,  used  in  radio  engineering, 
medicine  and  industrial  design,  have  received  proper 
recognition.   Fueled  portable  thermoelectric  generators, 
especially  for  military  applications,  are  already  used. 

In  the  present  work  an  attempt  is  made  to  design,  at 
least  theoretically,  a  thermoelectric  air  cooling  device 
servicing  the  passenger  compartment  of  a  small  automotive 
vehicle,  in  conjunction  with  a  thermoelectric  generator 
using  the  waste  heat  of  the  exhaust  gasses  of  the  engine 
to  provide  a  portion  of  the  required  electrical  power. 


The  feasibility  of  the  replacement  of  the  conventional 
alternator  by  a  thermoelectric  generator  is  also  examined. 

The  introduction  of  the  thermoelectric  generator  will 
improve  the  coefficient  of  performance  of  the  thermoelectric 
cooler,  and  will  also  increase  the  initial  installation  cost. 

The  design  was  developed  from  the  macroscopic  point  of 
view,  considering  the  present  state  of  the  art  of  thermo- 
electric materials. 

The  total  work  is  divided  into  five  sections: 

1)  The  heat  flow  analysis  for  the  semiconductor  material 
under  a  temperature  gradient  and  internally  generated  heat 
due  to  the  flow  of  electric  current.   The  results  of  this 
analysis  are  almost  common  to  the  thermoelectric  cooler  and 
generator.   In  this  section  a  differential  equation,  which 
describes  the  heat  balance  in  the  thermoelement  has  been 
derived.   In  this  equation,  the  cross  section  of  the  thermo- 
element is  treated  as  a  function  of  the  length  and  the 
properties  of  the  material,  as  thermal  conductivity,  Thomson 
coefficient,  etc.,  are  treated  as  functions  of  temperature. 

2)  The  analysis  and  design  of  the  thermoelectric  cooler. 
A  detailed  analysis  of  the  T.E.  cooler  is  made  for  both  modes 
of  operation,  maximum  heat  pumping  rate  and  maximum  coeffi- 
cient of  performance.   The  analysis  is  closed  with  comparison 
between  the  two  modes  and  the  selection  of  max  COP  as  the 
most  suitable  for  the  present  application. 

The  analysis  is  followed  by  an  indicated  design  of  a 
cooler  with  a  capacity  of  7000  BTU/hour  which,  from  a  quick 


estimation,  was  found  to  be  satisfactory  for  the  selected 
environmental  conditions . 

3)  The  analysis  and  design  of  the  thermoelectric 
generator.   The  analysis  is  conducted  from  the  standpoint 
of  maximum  efficiency  mode.   It  is  followed  by  selection  of 
the  proper  thermoelectric  materials  for  the  assumed  conditions 
Computer  techniques  are  employed  for  estimation  of  the 
temperature  relations  of  the  thermoelectric  material  proper- 
ties.  The  analysis  is  closed  with  the  estimation  of  the 
available  heat  power  from  the  exhaust  gasses.   The  attempted 
design  concerns  a  segmented  type  thermoelectric  generator 
which  is  carried  out  exclusively  by  a  computer  program. 

4)  Conclusions  and  recommendations  for  further  study. 

5)  Appendices 

Generally,  the  proposed  configuration  would  be  as  it  is 
shown  in  Figure  1. 
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FIGURE  I 

DIAGR/iT.lMATTC   ARR.ANGEI5ENT   OF  THE   PROPOSED   SYSTEM 
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II 


STEADY  STATE  HEAT  FLOW  IN  THE  SEMICONDUCTOR  ELEMENT 


The  derivation  of  the  heat  flow  balance  equation  through 
the  thermoelement  is  presented  when  the  cross  section  of  the 
element  is  a  function  of  the  length  and  the  characteristics 
of  the  material,  electrical  resistivity  and  thermal  conduc- 
tivity are  functions  of  temperature. 

Under  the  proposed  geometry  of  the  semiconductor,  the 
one  dimension  heat  flow  is  evident. 


FIGURE  2 
PROPOSED  GEOMETRY  OF  THERMOELECTRIC  SEGMENT 
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The  semiconductor  is  considered  to  be  adiabatically 

insulated  on  all  sides  except  side  ABCD  and  EFGH. 

The  temperature  of  the  solid  is  assumed  to  be  T„  at  the 

H 

front  side  (ABCD)  and  T  at  the  rear  side  (EFGH) ,  when  T 
is  greater  than  T^.   The  characteristics  of  the  material 
are  considered  to  be  variable  and  more  specifically  to  be 
functions  of  temperature  which  can  be  expressed  as: 


1)  Conductivity  Ji(T)  =  a^  +  a^^T  +  a2T^ 

2)  Resistivity  P  (T)    =  b^  +  b^^T  +  b2T^ 


3)   Seebeck  coefficient  o(  (T)  =  c^  +  c,T  +  C2 


t2 


4)   Thomson  coefficient  6"(T)  =  T  ^#^^  =  c^T  +   2c^^ 

Considering  the  part  of  the  semiconductor  which  is 
confined  between  the  cylindrical  surfaces  of  radii  r  and 
r+dr,  let  q(r)  and  q(r+dr)  be  the  rates  of  heat  flowing  in 
and  out  of  the  differential  element. 

If  a  current  I  is  flowing  through  the  semiconductor 
Joule  heat  will  be  generated  due  to  the  electrical  resistance 
of  the  semiconductor  and  heat  will  be  released  or  absorbed  — 
depending  on  the  direction  of  the  flowing  current  —  due  to 
the  Thomson  effect.   The  differential  of  the  Joule  heat 
generation  rate  is  given  as : 


q  =  t^^^r  (1-1) 

\      ^    A(r) 
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where  A(r)  is  the  surface  of  the  differential  solid  element 
at  radius  r. 

Also  the  differential  form  of  the  Thomson  effect  is 
given  as: 

^  ^  l^CDdT  (1-2) 

Consider  now,  for  balance  the  rate  of  heat  flow  out  of 
the  differential  solid  element  must  be  equal  to  the  rate  of 
heat  flow  in,  plus  the  rate  of  heat  generation  in  the  element. 
These  can  be  expressed  as: 

\W^^^)  =  '\(r)*'\^^^^  (1-3) 

But  q(r+dr)  can  be  expanded  in  a  series  of  the  form: 

If  the  first  two  terms  will  be  considered,  then  equation 
(1-3)  can  be  written  as: 

^lil)dr  =q  ±q  (1-5) 

Substituting  the  equivalent  relations  for  qj  and  q   equation 
(1-5)  becomes 

JM-^f_eciLiI^cT)^  (1-6) 
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The  conduction  heat  —  rate  of  flow  —  q(r)  is  given  by: 

Cj^(r):r  -  n^CT)>A(r)-^  (1-7) 

By  differentiation  in  respect  to  r  equation  (1-7)  becomes: 


\  (1-8) 


By  sxibstituting  of  equation    (1-8),    equation    (1-6)    becomes: 

:i(T)AM^4MA(r}[iljl^(T)dMJdLii^(T)^jma^ 


For  the  geometric  configuration  of  Figure  2  the  cross 

lection  area  A(r)  is  a  linear  function  of  radius  r: 

ACO  -  ^xr  (1-10) 

ACr)^^^  (l-lOa) 


where  m  =  ^'X  . 

Using  the  above  expression  for  the  cross  section, 
equation  (1-9)  becomes: 
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Equation  (1-11)  is  the  differential  equation  which 
describes  the  rate  of  heat  flow  through  the  semiconductor 
material,  with  the  geometry  of  Figure  2  and  the  conditions 
stated  before,  when  a  temperature  gradient  is  impressed  and 
cin  electrical  current  is  allowed  to  flow. 

This  equation  is  non-linear  and  at  least  a  rigorous 
analytic  solution  is  not  obvious. 

If  in  the  range  of  temperature,  when  the  thermoelectric 
couple  is  operated,  the  characteristics  of  the  semiconductor 
material  are  constant  or  reasonably  can  be  assumed  to  be 
constant,  then  equation  (1-11)  is  greatly  simplified  and 
takes  the  following  form: 

Since  the  temperature  gradient  under  which  the  thermo- 
electric element  of  a  cooler  is  operating  is  small,  equation 
(1-12)  reasonably  can  be  used  as  the  basic  equation  for  the 
design. 

This  is  not  the  case  for  the  thermoelectric  generator 
where  the  range  of  temperature  is  quite  extended,  and  the 
material  characteristics  present  sufficient  deviations  from 
the  corresponding  values  at  low  temperatures. 

For  this,  equation  (1-11)  will  be  used  as  the  basis  for 
the  thermoelectric  generator  design. 
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An  analytic  solution  of  equation  (1-11)  is  given  in 
APPENDIX  A  under  certain  assumptions.   APPENDIX  B  provides 
the  analytic  solution  of  equation  (1-12) . 
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III.   THERMOELECTRIC  COOLER 

A.   ANALYSIS  OF  THE  ELEMENTARY  COOLER 

The  operation  of  the  thermoelectric  cooler  is  based  on 
the  Peltier  effect,  which  was  discovered  in  the  eighteen- 
thirties.   This  phenomenon  was  explained  correctly  by  Lenz 
in  1838.   When  an  electric  current  I  flows  through  a  contact 
between  two  conductors,  a  certain  amount  of  heat  is  either 
evolved  or  absorbed  at  the  contact,  depending  on  the 
directions  of  the  current.   This  heat  is  known  as  the  Peltier 
heat  Q^. 

Becquerel  and  others  established  that  Qr,  ^^  directly 
proportional  to  the  current  I : 


Q.^  =  -nCT)I  (3-1) 


where  TTCt^  is  the  so-called  Peltier  coefficient. 

It  has  also  been  shown  that  if  the  current  flowing 
through  the  junction  has  the  same  direction  as  heat  flow 
which  is  due  to  the  heating  of  this  junction  then  the 
Peltier  heat  is  absorbed;  while  in  the  opposite  case,  heat 
is  evolved. 

The  Peltier  effect  has  been  observed  in  a  great  variety 
of  materials,  and  the  Peltier  coefficient  has  been  found  to 
be  dependent  on  temperature  in  a  complex  fashion. 
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The  thermocouple  is  considered  under  steady  state  condi- 
tions with  its  hot  junction  at  a  temperature  '!!„  and  the  cold 
junction  at  T  . 

Electric  current  is  considered  flowing  through  the  couple 
in  the  direction  opposite  to  its  Seebeck  voltage. 

The  arrangement  is  shown  in  Figure  3. 

The  law  of  conservation  of  energy  gives: 


V\-*  (3-2) 

where  q   is  the  rate  of  heat  absorption  at  the  cold  junction, 
q^  is  the  rate  of  heat  which  is  evolved  at  the  hot  junction, 
and  W  is  the  electric  power. 

An  important  efficiency  index  of  the  thermocouple  is  the 
coefficient  of  performance  which  is  defined  by: 


COP  =  -^  (3-3) 


Another  important  index  is  the  heat  pumping  rate  —  or 
the  refrigerating  effect  —  q  which  can  be  expressed  by: 


V^a)I  .^/^)^M^       ^XWA.(T)^ 


(3-4) 


The  thermoelectric  coolers  are  usually  operating  under 
a  small  temperature  difference.   This  is  particularly  true  in 
the  air  conditioning  applications.   Under  these  small 
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FIGURE   3 

ELEMENTARY   THERI,!OELECTRIC    COUPLE 
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temperature  ranges,  the  variations  of  the  properties  of  the 
semiconducting  materials  are  limited,  and  the  use  of  average 
values  introduces  slight  errors  of  the  order  of  one  percent. 

The  use  of  average  values  simplifies  the  analysis  of  the 
system  and  they  are  used  throughout  this  section. 

From  APPENDIX  B  the  rate  of  heat  flowing  to  the  cold 
junction,  because  of  the  existing  temperature  gradient  and 
the  generated  Joule  heat,  for  the  P-type  and  N-type  segments 


-'W^^'^ 


a  AC*-)  — 


r.Ci 


-«.i 


(3-5) 


(3-6) 


Introducing  the  following  variables 


)<  = 


C4t) 


(3-7) 


:}-- 


(3-8) 


n(Tc)=  Srav  (U/T,)UU,a)|)Tc 


(3-9) 
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equation  (3-4)  becomes: 

<^^-.  5(-t)i  T_  [^ .  ii]  AT  -  ^  r  [  f/  -  a]         (3-10) 

The  work  done  by  the  electric  source,  per  unit  time, 
consists  of  the  Joule  heat  evolved  in  the  thermocouple  arms, 
and  the  work  done  against  the  thermo-e.m.f  To(AT 

V/=  I^AT+I^H  (3-11) 

where  R  is  the  ohmic  resistance  of  the  couple.   The  coeffi- 
cient of  performance  can  be  written  as: 


COP-  n(Tc^J-[^'^T]AT-|rLf,>-^ej]  ^3.^,^ 


i^M  ^i^e.x-^U] 


Cop 


IRM  +  I^K 


(3-13) 


B.   CONDITIONS  FOR  MAXIMUM  HEAT  PUMPING  RATE 

The  heat  pumping  rate  of  the  elementary  thermoelectric 
cooler,  as  it  appears  in  equation  (3-10)  ,  is  a  function  of 
X,  Y,  and  I.   Relative  maximum  can  be  found  by  the  Lagrange 
multiplier  technique  using  as  constraint  function: 
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Thus  the  following  conditions  must  be  satisfied; 


9%o 

Af,. 

lO 

9% 

'^l- 

=  0 

(3-15) 


(3-16) 
9^o 
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=  0  (3-17) 


where  h  is  the  Lagrange  multiplier. 

Eliminating  h  from  equations  (3-15)  and  (3-16)  the  first 

condition  for  q   to  be  maximum  is  found  to  be: 
^o 


,     .  ^  .    .y. 


\    ^"^P  (3-18) 

By  substitution  of  equations  (3-7)  and  (3-8) ,  equation  (3-18) 
becomes : 


1  e.^p  Y 


(3-19) 


The  second  condition  for  the  maximum  heat  pumping  rate  is 
found  from  equation  (3-17) : 

T  4  -Us 

■*-°?^"  ^  (3-20) 

The  voltage  to  be  applied  to  the  thermocouple  is  given  by: 


V=SAT-tlo^t^  (3-21) 
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Substituting  the  expression  for  I  from  equation  (3-2)  in 
the  above  equation  it  is  found: 

The  required  electrical  power  for  the  couple  is : 

p=  VI  =  ^T^IoJ>^  .  (3-23) 


The  maximum  heat  pumping  rate  q    is  found  from  equation 

max 

(3-10)  after  the  substitution  of  equations  (3-18)  and  (3-20) . 
^m^^         2R.  (D-l)(D  +  l)  *^'^^' 


where 


'k 


^-     [^^%\  (3-25) 


___cCi____    T^-fTc 


(3-26) 


f/'f 


T_  Jh  (3-27) 
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The  coefficient  of  performance  under  the  maximum  heat  pumping 
rate  conditions  is  found  from  equation  (3-13)  after  trans- 
forming it  by  means  of  equations  (3-18)  and  (3-20) . 

COP       CO-3)CD^3)  (3-28) 

The  rate  of  heat  rejection  is  determined  with  the  help  of 
equations  (3-3),  (3-28)  and  (3-24). 


(3-29) 


C.   CONDITIONS  FOR  MAXIMUM  COEFFICIENT  OF  PERFORMANCE 

As  it  was  in  the  case  of  maximum  heat  pumping  rate,  the 

coefficient  of  performance  equation  (3-13) ,  is  also  a  function 

of  X,  y,  and  I. 

Applying  the  same  technique,  the  following  conditions 

must  be  satisfied: 


(3-30) 


9k 

9x 

..^=o 

9K 
9^ 

-+   \\P  -  o 

% 

=  o 

(3-31) 


(3-32) 


^I 
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Where  the  notation  K  instead  of  COP,  for  the  coefficient  of 
performance  is  introduced  for  convenience,  and  h  again 
represents  the  Lagrange  multiplier.   Elimination  of  h  between 
equations  (3-30)  and  (3-31)  gives  exactly  the  condition  of 
equations  (3-18)  and  (3-19) . 

By  algebraic  manipulation  the  following  expressions  are 
obtained: 


T-    °^  ^'^  (3-33) 

(D-l)K 

R&TD  (3_34) 

(D-l) 


V= 


Th 


(3-35) 


It  is  obvious  that  equation  (3-2  5)  goes  to  zero  when: 


D=—  (3-36) 


and  it  goes  to  infinity  when: 

AT  =  0  ^3-37) 

Combining  equations  (3-10),  (3-18)  and  (3-33),  the  heat 
pumping  rate  becomes: 


«^D .  (t,  -  T^)  (7^D  -  T«)       ( 3-38) 
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D.   COMPARISON  OF  THE  EXTERNAL  CONDITIONS 

The  two  cases  of  maximum  heat  pumping  rate  and  maximum 

coefficient  of  performance  differ  very  markedly  as  regards 

both  the  working  conditions  —  voltage  required  —  and  the 

efficiency  indices . 

Comparing  the  values  of  q    and  q   in  Figure  4,  it  is 

max      o 

seen  that  q„^^  is  greater  than  q^  and  consequently,  in  order 
max  o 

to  have  the  same  heat  pumping  rate  for  thermocouples  working 

under  conditions  of  q    and  COP    it  is  necessary  to  use  a 
max        max  -^ 

larger  number  of  thermoelectric  units  in  the  second  case  than 
in  the  first.   This  aspect  is  not  decisive  because  from  the 
efficiency  point  of  view  it  is  first  of  all  necessary  to 
ensure  that  the  coefficient  of  performance  is  sufficiently 
high,  since  it  is  this  coefficient  which  governs  the  elec- 
trical power  consumed  in  a  cooling  device.   On  the  other  hand, 
it  follows  from  Figure  5  that  the  coefficient  of  performance 

at  the  COP    is  greater  than  COP  and  that  this  difference 
max 

can  not  be  compensated  at  all.   The  preferable  conditions 
are  those  when  the  coefficient  of  performance  is  a  maximum 
since  the  materials  known  at  present  do  not  give  such  a  high 
coefficient  of  performance  that  it  can  be  afforded  to  have 
it  lower,  even  for  the  purpose  of  economy  of  material. 

Designs  based  on  the  maximum  heat  pumping  rate  could 
be  used  only  for  the  rapid  cooling  of  an  enclosure  or  for 
cooling  it  to  a  lower  temperature  at  the  cost  of  additional 
electrical  power  consumption. 
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E.   ESTIMATION  OF  TEMPERATURES  OF  OPERATION 

Complete  air  conditioning  involves  the  simultaneous 
control  of  the  following  factors: 

1)  Temperature 

2)  Humidity 

3)  Air  motion 

4)  Air  distribution 

5)  Dust 

6)  Bacteria 

7)  Odors 

8)  Toxic  gasses 

Of  these  the  first  three,  temperature,  humidity,  and 
air  motion,  are  paramount,  and  no  air  conditioning  system 
deserves  the  name  unless  it  satisfactorily  controls  those 
three  factors . 

In  the  present  work  the  thermoelectric  air  conditioning 
is  examined  only  from  the  standpoint  of  the  air  cooling 
process . 

The  device  should  be  also  considered  as  equally  appli- 
cable for  air  heating  process  by  just  reversing  the  direction 
of  the  current  through  the  thermoelements . 

From  the  standpoint  of  human  comfort,  the  object  of  the 
air  cooling  device  is  not  to  cool  the  body  but  to  control 
its  rate  of  heat  loss  within  reasonable  limits  by  methods 
that  are  known  to  be  the  most  conductive  to  comfort. 

There  are  four  factors  which,  by  their  combined  effect, 
determine  the  cooling  power  of  the  atmosphere  upon  the  human 
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body.   They  are  temperature,  humidity,  air  motion  and 
radiation. 

The  effect  of  humidity  and  air  movement  are  matters  of 
common  experience.   When  the  indoor  air  is  very  dry,  a 
higher  temperature  is  necessary  for  comfort  because  of  the 
greater  cooling  effect  of  the  dry  air.   The  evaporation  of 
the  perspiration  with  which  the  human  body  is  always  par- 
tially covered  has  a  marked  influence  upon  the  rate  of  heat 
loss,  and  the  rate  of  evaporation  is  dependent  upon  the 
h\imidity  of  the  surrounding  atmosphere. 

The  cooling  action  of  air  motion  upon  the  skin  increases 
heat  removal  by  the  increased  evaporative  effect  and  by  an 
increase  in  the  convection  effect  as  the  layer  of  air  in 
contact  with  the  skin  is  repeatedly  renewed. 

Loss  of  heat  by  radiation  from  the  exposed  body  and 
clothing  depends  on  the  temperature  of  the  surrounding  walls 
and  is  independent  of  air  temperature. 

Air  movement  reduces  to  some  extent  the  radiation  loss 
by  lowering  the  surface  temperature  of  the  body.   The  rate 
of  heat  loss  varies  according  to  Steffen's  fourth  power  law 
and  is  dependent  upon  the  difference  in  temperature  between 
the  surface  of  the  body  and  the  mean  temperature  of  the 
environment.   This  latter  temperature  is  referred  to  as  the 
mean  radiant  temperature. 

Experiments  made  on  a  large  number  of  human  subjects 
have  demonstrated  that  there  is  a  consistent  relation  between 
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the  temperature,  humidity  and  air  motion  in  their  effect 
upon  comfort.   When  the  results  of  such  observations  are 
plotted  on  a  psychometric  chart,  it  is  found  that  the  points 
representing  conditions  of  equal  warmth  fall  in  approximately 
straight  lines.   These  straight  lines  are  called  "effective 
temperature  lines".   For  moving  air  the  effective  temperature 
for  a  given  temperature  and  humidity  is  lower  than  that  corres- 
ponding to  still  air.   The  effective  temperature  is  the 
only  true  index  of  human  comfort  that  is  now  available. 
Neither  the  dry-bulb  nor  the  wet-bulb  thermometer  alone 
trully  indicates  the  conditions  of  the  air  from  a  comfort 
standpoint. 

The  optimum  conditions  for  human  comfort  at  summer  time 
are  considered  to  be  71  F°  effective  temperature  and  air 
movement  of  15  to  35  feet  per  minute  for  an  outdoor  temperature 
of  80  F**. 

The  following  table  gives  the  corresponding  effective 
temperatures  for  several  outdoor  dry-bulb  temperatures. 


OUTSIDE  DRY-BULB  EFFECTIVE 

TEMPERATURE  (F°)                TEMPERATURE  (F°) 

105  75.5 

100  75 

95  74 

90  73 

85  72 

80  71 


TABLE  1 
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These  effective  temperatures  are  the  statistical  results 
of  extensive  experimental  work,  where  several  factors  were 
considered,  between  them  the  stress  which  a  human  body  suffer; 
when  it  is  moving  from  low  to  high  temperature  environments 
euid  vice  versa. 

The  outdoor  conditions  to  be  assumed  in  the  calculations 
of  the  cooling  and  dehumidifying  load  should  not  be  the 
highest  recorded  temperature  in  a  location  which  usually 
exist  for  only  a  few  hours  during  the  year.   They  should  be 
chosen  as  conditions  which  may  be  exceeded  for  a  few  days 
during  the  summer. 

In  every  large  city  a  certain  design  dry-bulb  and  wet- 
bulb  temperature  which  are  readily  ascertained,  have  become 
established  and  generally  adopted.   The  following  table  gives 
these  temperatures  for  a  few  cities  of  the  U.S.A. 


LOCATION 

New  York  City 

Chicago 

Philadelphia 

Detroit 

Cleveland 

St.  Louis 

Washington 

New  Orleans 

San  Francisco 

Dallas 

Jacksonville 


DESIGN  TEMPERATURES 


DRY -BULB 

WET-BULB 

95 

75 

95 

75 

95 

78 

95 

75 

95 

75 

95 

78 

95 

78 

95 

79 

90 

65 

100 

78 

95 

78 

TABLE  2 
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since  an  automobile  is  not  designed  to  operate  in  a 
certain  place,  the  worst  case  outdoor  conditions  are  selected 
as  those  corresponding  to  the  city  of  Dallas. 

So,  outdoor  dry-bulb  temperature  is  100  F°  and  indoor 
maintained-effective-temperature  the  corresponding  75  F**. 
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IV.   COOLER  DESIGN 

A.   GENERAL  CONSIDERATIONS 

The  design  is  based  on  the  worst  case  which  corresponds 
to  an  outdoor  temperature  of  100  F**  and  an  effective  tem- 
perature, for  human  comfort,  of  75  F°.   This  means  that  the 
temperature  of  the  incoming  air  in  the  hot  side  of  the 
thermo  elements  would  be  100  F**  and  the  maintained  tempera- 
ture in  the  passenger's  compartment  of  the  vehicle  would 
be  75  F". 

The  proposed  configuration  of  the  thermoelectric  cooler 
is  shown  in  Figures  6a,  6b.   It  consists  of  radially 
arranged  thermocouples,  forming  rings.   A  number  of  such 
rings  are  assembled  to  form  the  annular  configuration.   Fins 
have  been  introduced  in  both  sides  of  the  thermocouples  to 
increase  the  rate  of  heat  transfer.   The  outer  annular  finned 
space  forms  the  heat  absorbing  exchanger  and  the  inner  cylin- 
drical finned  space  forms  the  heat  rejecting  exchanger.   The 
cooling  medium  is  air  which  flows  by  means  of  fans,  one  for 
every  heat  exchanger.   When  the  vehicle  is  moving  with  a 
speed  higher  than  a  certain  level,  there  is  no  need  for 
operation  of  the  heat  rejecting  exchanger  fan. 

It  is  also  assumed  that  the  vehicle  has  been  parked  for 
a  long  period  under  the  sun  and  the  temperature  of  the  indoor 
air  has  reached  100  F° , 

The  temperatures  of  the  hot  and  cold  side  of  the  thermo- 
elements, T  ,  T  ,  can  be  estimated  if  the  desired  logarithmic 
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FIGURE  6a 
CROSS  SECTION  OF  T.E.  COOLER 
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FIGURE  6b 
LONGITUDINAL  CROSS  SECTION  OF  T.E.  COOLER 
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mean  and  output  temperatures  of  the  cooling  air,  in  both 
sides  of  the  thermocouples  are  selected. 

According  to  the  definition  of  mean  logarithmic 
temperature: 


A"T  -    ^triav  "  ctm.iV  (4-1) 


the  hot  and  cold  shoe  temperatures  of  the  thermocouples  are 
estimated  as: 


Tcrt^       AT,^l-T,-^  (4-2) 

A-l 


where; 


AT 

=  e  K 


At      =  t.     -  t     ^ 
in  out 


desired  value  of  logarithmic  mean  temperature 


fi  .  ,  9     =   temperature  differences  as  shown  in 
^min  ^max     ^.^j;^^  ^^ 


The  cooling  capacity  of  the  thermoelectric  cooler  is 
selected  to  be  7000  BTU  per  hour.   This  figure  is  considered 
as  a  moderate  value,  sufficient  to  serve  a  small  vehicle 
carrying  four  passengers  under  the  climatological  conditions 
described  previously. 
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FIGURE  7 

ASSUMED  TEMPERATURE  DISTRIBUTION  FOR 
DETERMINATION  OP  LOG.  MEAN  TEMPER/.TURES 
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As  the  most  suitable  thermoelectric  material  for  cooling 
application  selected  a  bismuth  —  tellurium  —  antimony  com- 
pound for  the  P-type  semiconductor,  and  a  bismuth  —  tellutium 
selenium  compound  for  the  N-type.   These  materials  have  been 
developed  since  196  8  from  Westinghouse.   The  average  proper- 
ties of  the  materials  are  displayed  in  Table  3. 


ITEM 

UNITS 

P-type 

N-type 

Seebeck  Coeff. 

V/deg  K 

205  X  lO"^ 

-  220  X  lO"^ 

Resistivity 

ohm- cm 

1.2  X  10"^ 

1.08  X  lO"^ 

Th.  Conductiv. 

watts/deg 

c-cm 

0.0121 

0.0149 

Fig.  of  merit 

1/deg  K 

2.82  X  lO'^ 

3.0   X  lO"^ 

TABLE  3 

In  order  to  obtain  good  mechanical  strength  the  materials 
are  processed  by  pressure  and  sintering. 

B.   COOLER  DESIGN  BASED  ON  MAXIMUM  COP 

This  design  is  based  on  the  results  obtained  from  the 
analysis  in  Section  III-C. 

1.   Thermoelectric  Couple 

The  coefficient  of  performance,  as  was  defined  in 
Section  III-C,  is  independent  of  the  geometric  configuration 
of  the  thermoelectric  couple.   It  is  only  a  function  of  hot 
and  cold  shoe  temperature  and  of  the  material  properties. 
The  heat  pumping  rate  and  the  optimum  current  depend  on  the 
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total  electrical  resistance  of  the  thermocouple  which  is 
a  strong  function  of  geometry.   From  this  standpoint,  an 
optimization  of  the  geometry  had  been  attempted  using 
iterative  computer  techniques  and  the  selected  dimensions 
in  Table  4. 


ITEM 

UNITS 

P-type  segment 

N-type  segment 

^1 

cm  ( in . ) 

7.0  (2.756) 

7.0  (2.756) 

R2 

cm  (in.) 

8.5  (3.346) 

8.5  (3.346) 

<^ 

degrees 

10.0 

11.0 

X 

cm  (in.) 

1.0  (0.3937) 
TABLE  4 

1.0  (0.3937) 

The  calculated  cooling  and  electrical  characteristics 
of  the  thermocouple  are  shown  in  Table  5. 


ITEM 
Electrical  Resistance 
Hot  shoe  temp.  T„ 
Cold  shoe  temp.  T^ 
Electrical  current 
Applied  Voltage 
COP 

Required  Power 
Heat  pumping  rate 
Heat  rejection  rate 


UNIT 
ohm 
degr.  K 
degr.  K 
amperes 
volts 

dimension less 
watts 
watts 
watts 


VALUE 

2.31  X  10 

322.73 

291.5 
14.4 
0.048 
1.13534 
0.691 
0.74466 
1.43569 


-3 


TABLE  5 
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For  the  cooling  load  of  7000  BTU  per  hour  a  total 
number  of  2754  thermocouples  is  required.   If  a  nine  degree 
wedge  of  insulation  is  provided  for  each  thermocouple, 
then  12  couples  form  one  ring  and  the  total  number  of  rings 
is  determined  as  230.   Each  ring,  according  to  Table  4, 
has  width  1.0  cm.   Allowing  insulation  of  0.3  cm  between 
the  rings,  the  total  length  of  the  configuration  is  299  cm. 

Four  units,  each  having  a  length  of  74.75  cm  and 
consisting  of  57.5  rings,  will  satisfy  the  requirements. 
The  externally  supplied  electrical  power  is  calculated 
to  be: 

2754  X  0.691  =  1903  watts 

and  the  heat  rejection  rate: 

2754  X  1.43569  =  3954  watts  or  13504  BTU/hour 

2.   Cold  Side  Heat  Exchanger 

The  introduction  of  extended  surfaces,  in  the  form 
of  fins  with  rectangular  cross  section,  increases  the 
exchange  surface  of  the  heat  exchanger.   The  fins  have  an 
efficiency  given  by: 


^        -_         -IcxnrvlvCmlD^  (4-3) 
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where 


om. 


Bit 


(4-4) 


Figure  8b  shows  the  variation  of  n^  as  a  function  of 
for  several  values  of  h.   The  cooling  medium  is  air  with 
average  properties  displayed  in  Table  6. 


SYMBOL 


ITEM  UNITS  VALUE 

Thermal  conductivity  BTU/f t-hour-F"   0.015 

Dynamic  viscosity  Ibm/ft-sec.      1.7  x  10 

Density  Ibm/ft^         0.075 

Thermal  capacity  BTU/lbm-F*'      0.24 


-5 


TABLE  6 

Figure  9  shows  the  heat  exchanger  arrangement.   The  long 
fins  form  four  different  types  of  air  ducts,  namely  A,  B, 
C,  and  D. 

Each  air  duct  has  an  equivalent  hydraulic  diameter 
given  by : 


e    p 


(4-5) 


where  S  is  the  cross  section  and  P  the  perimeter  of  the 
duct. 
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FIGURE  9 
ELEMENTARY  SECTION  OF  COLD  SIDE  HEAT  EXCHANGER 
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The  calculated  equivalent  diameters  for  each  duct 
type  are: 


(D^)^  =  48.284  X  10~^  ft 


(D^)g  =  43.89   X  lO"^  ft 


(Dg)^  =  43.89   X  lO""^  ft 


(Dg)j^  =  62.26   X  lO"^  ft 


The  selected  temperatures  for  the  air  and  thermo- 
couple cold  shoe  are  displayed  in  Table  7. 

ITEM  TEMPERATURE 

Incoming  hot  air  100  F° 

Outgoing  cold  air  70  F° 

Maintained  room  environment  75  F° 

Thermocouple  cold  shoe  291.5  K** 

TABLE  7 

For  operation  at  low  noise  level,  air  velocity  is  selected 
to  be  1800  feet  per  minute. 

The  Reynold's  number  is  determined  from: 


NR^  =  -LDeJVL_  (4.6, 
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Using  the  air  properties  of  Table  7  the  corresponding 
value  of  N   for  each  type  of  duct  is  calculated  as : 


<'^re'a=  32S5 


(Nre'b  =   8""° 


(Nre'c  =   8"" 


(Nre'd  =  1"" 


These  values  are  indicative  of  turbulent  flow.   The  calcu- 
lation of  the  convection  heat  transfer  coefficients  is 
based  on  the  Dittus  and  Boelter  equation: 


N^„=  0.0265  (Nj^)°-8  (V'°'^  '"■" 


where  N....  and  N   are  the  Nusselt's  and  Prandtl's  number 
NU      pr 

correspondingly,  given  by: 


W   -  ^^g-  (4-8) 


^r=^  (4-9) 

For  a  Prandtl's  nuinber  almost  equal  to  0 . 7  the  Nusselt's 
number,  for  the  different  type  of  air  ducts,  is: 
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("nu'a 

- 

35.1428 

<«nu'b 

= 

32.56 

%u'c 

= 

32.56 

%u'd 

= 

43.07 

The  corresponding  convection  heat  transfer  coefficients  are 
found  from  Equation  (4-8)  as: 


*>A 

= 

10.92 

hfi 

= 

11.128 

he 

= 

11.128 

hj^  =   10.376 


Using  the  above  values  for  the  convection  heat  transfer 
coefficient  and  assuming  aluminum  as  the  fins  material,  with 
thermal  conductivity  92  BTU/f  t-hour-F** ,  the  efficiency  of 
the  long  and  short  fins,  according  to  Equation  (4-3),  is: 


n^,   =   0.9256    and   n^^   =   0.9716 


correspondingly . 

The  rate  of  heat  absorption  is  calculated  from 
the  following  equation: 


49 


0,=    Ml^T.V:^S^  (4-10) 

where  the  subscript  j  corresponds  to  different  types  of  air 
ducts  of  the  heat  exchanger.   Analytically  the  exchange 
effective  surfaces  per  thermocouple  column  are: 


^A  "  ^^"fl  "^  "^"fs  "^  -^'^^^    ^^°  "  1983.2  cm^  =  2.135  ft^ 

Sg  =  (5n^^  +  0.26)  230  =  1124.24  cm^  =  1.21  ft^ 

^C  "  5n^j^230  =  1064.44  cm^  =  1.146  ft^ 

^D  "  ^^^fl  "^  -^^fs  "^  1-232)  230  =  2018.2  cm^  =  2.172  ft^ 


From  the  temperatures  of  Table  2  the  logarithmic  mean 
temperature  of  the  cooling  air,  according  to  Equation  (4-1) , 
is  found  to  be: 


AT(^=  15.0 


Substitution  of  the  determined  values  for  S .  and  h.  in 
Equation  (4-10)  gives: 


q  =  (2.135-10.92  +  1.21-11.128  +  1.146-11.128 

+  2.172-10.376)  A"E,   =  1081  BTU/hour 
In 
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For  the  12  thermocouple  columns  of  the  device,  the  capacity 
of  the  heat  absorbing  exchanger  is  found: 

q^Q^  =  12  792  BTU/hour 

The  above  figure  for  the  heat  exchanger's  capacity  offers 
a  margin  factor  of  1,85  which  ensures  satisfactory  operation 
under  even,  not  considered,  worse  conditions. 

The  air  flow  requirements  can  be  determined  from: 


m  = rrr^ r r-  (4-11) 

^p  ^^in  -  ^out^ 


for  q  =  12,792  BTU/hour,   c   =0.28  BTU/lbm-F° , 

t.^=   100  F° ,    and  t   ^  =  70  F*',  the  rate  of  air  flow  is 

determined  as: 


m  =  152  2.8  Ibm/hour 


For  the  air  density   f  =  0.075  Ibm/ft^ 


m  =  5.64  ft  /sec 


3.   Hot  Side  Heat  Exchanger 

This  is  also  a  finned  type  heat  exchanger.  It  con- 
sists of  a  tube  with  radially  introduced  fins  of  two  differ- 
ent lengths.   The  fin's  material  is  assumed  to  be  aluminum 
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with  thermal  conductivity  92  BTU/hour-f t-F° .   Following  the 
same  procedure  as  in  the  cold  side  heat  exchanger  the  obtained 
results  are  displayed  in  Table  8. 


ITEM 

UNIT 

VALUE 

Required  capacity 

BTU/hour 

13,503 

Hot  shoe 

po 

121.5 

Incoming  air 
Temperature 

Outgoing  air 

po 
po 

100 
110 

Log.  mean  difference 

po 

16 

Tube  diameter 

ft 

0.4462 

Long  fin  length 

ft 

0.1833 

Short  fin  length 

ft 

0.082 

Niomber  of  air  ducts 

- 

24 

Duct  equivalent  diameter 

ft 

0.03784 

Duct  length 

ft 

9.55 

Air  velocity 

ft/sec 

86 

Reynold ' s  number 

- 

20,860 

Nusselt's  number 

- 

67.15 

Prandtl's  number 

- 

0.7 

Convection  coefficient 

BTU/ft^-hour-F** 

26.6 

Long  fin  efficiency 

- 

0.6296 

Short  fin  efficiency 

- 

0.882 

Calculated  capacity 

BTUAour 

18,983 

Security  factor 

- 

1.4 

Volumetric  rate  of  flow 

ft^/sec 

9.4 

TABLE  8 
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FIGURE  10 

ELEWENTARY  T.C.-HOT  SIDE 
HEAT  EXGUAMGER 
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V.   THERMOELECTRIC  GENERATOR 

A.   ANALYSIS  OF  ELEMENTARY  GENERATOR 

The  elementary  thermoelectric  generator,  which  is  operating 
under  steady  state  conditions,  consists  of  a  thermocouple 
externally  connected  to  a  resistive  load,  as  it  is  shown  in 
Figure  11. 

The  operation  of  the  device  can  be  described  as  follows: 
the  hot  junction  of  the  thermocouple  is  assumed  at  temperature 
Tp.   The  thermocouple  arms  are  assumed  to  be  adiabatically 
insulated,  and  heat  exchange  takes  place  only  at  the  junctions. 
At  the  hot  junction  heat  q   is  absorbed,  and  heat  q,  is 
evolved  at  the  cold  junction  per  unit  time.   Under  these 
conditions,  a  thermo-e.m.f .  E  =  o^  AT  is  established  in  the 
circuit  and  it  maintains  an  electrical  current  I.   The 
coefficient  in  the  thermo-e.m.f.  is  the  Seebeck  coefficient 
of  the  couple 


<X  -    t^pCT) 


^(t) 


(5-1) 


From  the  law  of  energy  conservation  it  follows: 


\-    \^   ^  '^-^' 


where  W  is  the  available  electrical  power  at  the  external 
load. 


54 


FIGURE  II 
ELEMENTARY  T.E.  GENERATOR 
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(5-3) 


where  R  is  the  total  electrical  resistance  of  the  thermocouple 

Equations  (5-3)  suggest  the  electrical  equivalent  of  the 
thermocouple  shown  in  Figure  12 . 


Aaaaa 


Figure  12:   Electrical  equivalent  of  thermoelectric 
couple. 


The  general  expression  for  the  heat  flov;  rate  q^  has  the 
following  form: 


5-4) 
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This  expression  is  reduced  to  a  more  practical  form  by 
introducing  the  estimate  of  the  rate  of  heat  flow  in  the 
conducting  solid,  as  it  is  determined  in  Appendix  A 


-.-^i^AT-fiiM^±iaT    ( 


5-5) 


where 


•c  -' 

.   j_  /6a)(/^;)CT)cJT)oiT 


(5-6) 


The  plus  or  minus  sign  in  equation  (5-5)  allows  for  two 
possible  directions  of  the  electrical  current:   the  minus 
sign  representing  the  current  in  the  P-type  segment  and  the 
plus  sign  giving  the  current  in  the  N-type  segment.   As  it 
is  shown  in  Appendix  A,  the  value  of  P  is  very  close  to  half 
of  the  average  resistivity  for  the  range  of  temperature 

For  convenience  in  further  calculations,  it  is  desirable 
to  replace  equation  (5-5)  with  a  less  exact  but  much  simpler 
expression: 
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cJr 


(ill/^T-lff-MlijaT   (5- 


where: 


C4x) 


7) 


Combining  equations  (5-4)  and  (5-7) ,  the  rate  of  heat 
flow  at  the  hot  shoe  of  the  thermocouple  becomes : 


1^  ^   >v  V  .  _     (5-8) 
(6  -6^)  AT 


Introducing  the  Kelvin's  relation,  equation  (5-8)  becomes: 


<^^,-,Ta.[-^.-^]AT-if[f,x-e.^j]-iLvMAT 


(5-9) 


where 


and 


P^ 


X  - 


Mi) 


(5-10) 


^--■^i^T(<'rtT))-(<<,CT))]olT 
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The  efficiency  of  the  elementary  thermoelectric  generator 
is  defined  as: 


«T„I  .  [3f  ^  ^  j  AT-  ll^[PpX  +  fj]  -1  (syi)  AT 


\--  _,  ,1  .,   ^:::,:"^i.,  _,. ..  ^     (^-n. 


If  I  and  At  are  considered  as  constants,  then  the  efficiency 
'^  is  a  function  of  x  and  y  which  determine  the  geometry  of 
the  P  and  N  type  segments. 

Figure  13  shows  the  three  dimensional  representation  of 
efficiency  as  it  is  interpreted  in  (5-11)  ,  for  a  certain 
temperature  range  and  S.G   as  thermoelectric  material. 

B.   CONDITIONS  FOR  MAXIMUM  EFFICIENCY 

Equation  (5-11)  is  a  function  of  geometric  variables 
X,  y  and  of  the  electrical  current  I,  for  a  given  temperature 
range.   Using  the  Lagrange  multiplier  technique  again,  the 
necessary  conditions  for  maximum  efficiency  are  defined  as: 


2T^ 


W.O    ^    O  (5-12a) 


^  ^Up  =  0  '   (5-12b) 

i2L   =0  (5-12C) 

where  h  is  again  the  Lagrange  multiplier  and  R  = 
is  the  constraint  function. 
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FIGURE  13 
THREE  DIMENSION  REPRESENTATION  OF  EFFICIENCY 
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Elimination  of  h  from  equations  (5-12a)  and  (5-12b)  gives 
the  already  foiind  condition 


^    n.f.^^/^ 


Equation  (5-11)  can  be  written  as: 


(5-13) 


\- 


^lAT-I^'R 


°<THl4[-2r^^^/^T-ll^R.-IC4-^^)AT 


(5-14) 


Multiplying  and  dividing  by  R  and  introducing  the  condition 
of  (5-13),  the  efficiency  equation  becomes: 


where 


(5-15) 


(5-16) 


Optimization  of  (5-15)  with  respect  to  X   gives  the 
condition: 


i 


^  AT 


(5-17) 
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where : 


(5-18) 


The  optimijm  current  is: 


-r    ,  _   S  AT 


The  optimum  effective  voltage  is : 

V=  °^^  ^"^  (5-20) 

When  the  conditions  (5-13)  and  (5-17)  are  satisfied,  the 
efficiency  reaches  its  maximum  value,  which  is  found  from 
equations  (5-15)  and  (5-17) : 


-n    =    il    _Uzl) ,  (5-21) 

^'"'''   Th  (q  +  i  \ 

The  electrical  power  dissipated  in  the  external  circuit  is 
found  from  equation  (5-19)  and  (5-20) . 


The  formulae  developed  above  show  that  the  efficiency  of 
a  thermocouple  and  its  power  output  increase  with  increase  of 
the  working  range  of  temperatures.   The  thermocouple  efficiency. 
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like  the  cooler  CO. P.,  depends  only  the  the  ratio  T^/T„ 

C   n 

and  the  quality  factor  G,  which  is  slightly  different  from 
the  quality  factor  D  for  the  thermoelectric  cooler. 

For  the  present  application  the  working  temperatures  of 
the  thermoelectric  generator  are  varying  as  the  temperature 
of  the  exhaust  gasses  varies  for  different  power  level 
operation  of  the  automobile  engine.   In  the  present  state  of 
material  development,  no  single  thermoelectric  material  is 
satisfactory  for  use  over  the  complete  temperature  range. 
Therefore,  it  is  necessary  to  use  different  materials  in  each 
temperature  range.   One  way  to  accomplish  this,  is  to  form 
each  leg  of  the  couple  by  joining  different  materials  in 
series;  that  is,  forming  a  sandwich  structure  of  the  appropriate 
thermoelectric  materials  for  each  temperature  range. 

C.   OPTIMIZATION  OF  THE  SEGMENTED  T.E.  GENERATOR 

In  Figure  14  two  P-type  materials,  P^^  and  ^2'    ^^^   connected 
in  series  to  form  a  positive  thermoelectric  arm  and  two  N-type 
materials,  N,  and  N2 /  form  the  negative  arm  of  the  thermocouple. 

The  optimization  problem  is  to  maximize  the  thermal 
efficiency,  subject  to  specification  of  the  interface  tempera- 
tures T   and  T   over  the  equivalent  length  of  thermocouple, 
X      y 

L  =  In  ( -^  )  ,   equivalent  cross  sectional  area  Ap=^(4l)  , 
and  the  temperature  averaged  material  properties. 

The  reciprocal  of  the  thermal  efficiency,  as  defined  in 
equation  (5-11),  is: 

i,    ^THl-^rTH-TOkpi-^(:TH-T^)K^,- j;f(gH-^R^,M(gp-6w)AT  (5.23) 
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<> 


FIGURE  14 
ELEMENTARY  SEGMENTED  GENERATOR 
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Introducing  the  concept  of  equivalent  length  and  cross  section 
for  the  P-type  and  N-type  segment  as: 


Ap=  (4>t)p 


(5-24) 


The  terms  in  the  numerator  of  equation  (5-2  3)  have  the 
following  interpretations : 


Kp,= 

-^'■t 

Kna  ■=. 

Rpi  = 

^f^  Af 

fl«l-. 

(5-25) 
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"'"at 


(5-25) 


The  electrical  current  I  in  equation  (5-2  3)  is  given  by: 


I  =  5-^^  (5-26) 

o 


where:   R  =  Rpi  +  Rp2  "^  ^1  "*"  ^2  =  ^  "*"  ^  is  the  total 
electrical  resistance  of  the  thermocouple  and  E  is  the 
circuit  voltage  pro(duced  by  the  cumulative  Seebeck  effect 
in  the  positive  and  negative  elements.   The  segmented  elements 
may  be  considered  as  batteries  connected  in  series. 
The  voltage  E  is  given  by: 

E  =  ^'  (Tjj  -  T^)  (5-27) 


where; 


^  =  -^  I  /  '<r,(^^dT  .jd,^^  ct)jt  y^"^,  (T)dT  ycK^^(T)ciT  j 
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Substituting  equation  (5-26)  into  (5-23),  the  latter  becomes: 


±  -  (TH-T>)kpi-^(TH-Tv^i<fA/i-^^">"H  R-egp  "  %(SL^u^Y-^      ~  C^ri'^^^^   R-*i^o   (5-28) 


With  the  variables  9^  and  9^  defined  by  the  equations: 

(5-29) 


and  multiplying  the  numerator  and  denominator  of  equation 
(5-28)  by  (R  +  R  )   ,   equation  (5-28)  becomes: 

V  E=»Uo  (5-30) 

« 
With  the  quantities  R^,  ^,    K^^,   K^^  defined  by  the 
equations: 


(5-31) 


<- 

R^Ap 

Rn  = 

RnA~ 

k'n 
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equation  (5-30)  becomes: 


(5-32) 


or; 


(5-33) 
Assuming  _^  _  £.  [  ^^]      ,      minimization  of  equation  (5-33) 

A 

with  respect  to   (J^\       ,    gives: 


li^L  =  0.  C  r;  -  9r  kU:  (A.)\  o 
30  '" 


(5-34) 
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From  equation  (5-34)  the  condition  is  obtained: 


«.i_8£j4^y/^  ,5.35, 


Applying  condition  (5-35)  to  equation  (5-33) ,  the  latter 
becomes : 

With  the  figure  of  merit  defined  as: 


(5-36) 


M.  (^rT. 


(5-37) 


4[(e,R;K;,)\(©.Cx-V'^] 

equation  (5-36)  becomes: 


1       ^MCt^.-tJ        «^'(-t,-t,) 


With  the  partial  derivative  of  (— )  with  respect  to  -^      set 
equal  to  zero,  the  following  condition  is  obtained: 


(5-38) 
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(5-39) 


Substitution  of  (5-39)  into  equation  (5-38)  gives  the  minimum 


or  the  maximum  efficiency. 


Summarizing,  the  maximization  problem  of  the  segmented 
generator  efficiency  can  be  stated  as  follows:   Given  the 
values  of  ^^j,  f  ij '  ^"^ij'   "^H'   "^C  \'    ^^  ^  determine 
the  values  of  A^,  t'-i'  ^'    ^o'  ^'  ^^^  ^  '   where   i  =  P,  N 
and  j  =  1,  2 . 

Solution  of  the  problem  requires  that  the  following 
equations  be  satisfied: 


j_    o<  CTm-Tc)  (5-40) 


(5-41) 


'  2. 

h  =         ("^^  "'"« (5-42) 
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Jk  -[Tm(^-^)    Y'z  (5-43) 
MM        o<'     J 

L=  Cpi  +  £p2=L*Cz  (5-44) 

1^  R-PI  +  R.pi+R„it  H^j  (5-45) 


Equations  (5-46)  and  (5-47)  represent  the  energy  balance  for 
the  interfaces  between  elements  P,,  P2  and  N,,  N2 .   They 
state  that  the  amount  of  energy  flowing  to  interface,  plus 
or  minus  the  amounts  of  energy  generated  or  absorbed  at  the 
interface  is  equal  to  the  amount  of  energy  flowing  from  the 
interface. 

D.   T.E.  GENERATOR  DESIGN 

The  design  of  the  segmented  thermoelectric  generator  is 
implemented  by  computer  aided  iteration  methods.   Such  methods 
have  been  coded  for  the  IBM-360  as  it  is  displayed  in  program 
Number  one.   The  procedure  consists  of  several  steps: 
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1)   For  the  selected  materials,  calculation  of  o^ . 


^i-i'^Li      ^°^  ^  "  ^'^  ^^^  ^  "  ^'^' 


'3  ■  fi: 


2)  Calculation  of  o(   9^  0^, 

3)  Initial  estimation  of  c«   ?    9        9  from: 


4)  Calculation  of  R^ ,  P^,  Kp^^,  Kp2  • 

5)  Calculation  of  -r— ,  from  Equation  (5-41)  and  for 
selected  value  of  A^  determination  of  A^. 

6)  Determination  of  figure  of  merit  M  from  Equation 

(5-37) . 

R 

7)  Calculation  of  -^  from  Equation  (5-42) . 

8)  Calculation  of  R  from  equation: 


^-  Ce.,  Ks  -  e  J.o  ^^  ^  (^U  ^  eJ-0  i^ 


9)   Calculation  of  current  I  from  Equation  (5-40) . 
10)   Calculation  of  both  sides  of  Equations  (5-46)  and 
(5-47) .   If  their  differences  are  greater  than  a 
predetermined  figure,  then  a  new  set  of  Ip/ /  X^t    ij 
estimated  and  the  procedure  is  repeated  until  the 
wanted  accuracy  is  obtained. 
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11)  Calculation  of  R^,,  R^,. 

12)  Calculation  of  efficiency  from  Equation  (5-38) . 

The  procedure  is  based  on  an  "in  advance"  selection  of 
the  interface  temperatures  T„  and  T  with  the  criterion  of 
safe  behavior  of  the  materials  at  the  extreme  operating 
temperatures . 

When  the  lengths  of  the  P  and  N  type  segments  are  cal- 
culated for  the  extreme  temperatures ,  the  problem  of  re- 
estimation  of  the  interface  temperatures,  for  different  hot 
junction  temperatures,  is  encountered  by  modification  of 
the  iteration  procedure  in  step  10.   That  is,  for  the  calcu- 
lated lengths,  T^  and  T^  are  varying  until  Equations  (5-46) 
and  (5-47)  are  satisfied. 

The  adverse  operating  conditions  of  the  thermoelectric 
generator  —  high  hot  junction  temperatures  and  extended 
operating  range  —  impose  restrictions  in  the  selection  of  the 
best  available  thermoelectric  material  for  energy  conversion. 

All  thermoelectric  materials  represent  the  best  of 
their  characteristics  in  a  quite  narrow  range  of  temperatures 
This  problem  is  partially  encountered  by  the  segmented  type 
generator  design,  where  the  one  material  acts  as  the  comple- 
ment of  the  other  in  regions  of  poor  characteristics. 

For  the  present  application,  the  selected  thermoelectric 
materials  and  some  of  their  properties  are  displayed  in  Table 
9. 
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ITEM  MATERIAL  MAX  OPERAT.  TEMPERATURE 


Ag  Sb  Te^  900  (K°) 


^1 

P2  25%  Bi2Te3  +  75%  Sb2Te3  600  (K**) 

Nj_  Pb  Te  (25%  Sn  Te)  900  (K°) 

N2  75%  Bi^'Te^   +   25%  Bi2Se  700  (K*') 


TABLE  9 

With  proper  dopants  the  maximum  operating  temperature  of  P, 
and  N,  materials  can  be  considerably  extended. 

The  thermoelectric  properties  of  the  materials  are 
displayed  in  Figures  15  through  21.   The  data  have  been 
taken  from  Reference  [1]  and  [2]  and  the  corresponding 
temperature  relations  developed  using  the  "best  fit  poly- 
nomial" technique  as  it  is  shown  in  Program  2. 

With  these  data  supplied,  the  computer  runs  through  the 
iteration  procedure  described  earlier,  and  solves  for  the 
optimiom  thermal  efficiency  and  the  optimum  design  parameters. 
The  optimum  solution  is  displayed  in  Table  10. 

The  maximum  and  minimum  values  for  efficiency,  current, 
and  voltage  in  Table  10  are  associated  with  the  highest  and 
lowest  hot  junction  temperature. 

Figures  22  through  2  4  shows  the  variation  of  efficiency, 
voltage,  and  current  with  respect  to  hot  junction  temperature 
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ITEM 

UNIT 

VALUE 

REMARKS 

^H 

K° 

720 

min  -  940  max 

Selected 

''C     a 

K» 

335 

Selected 

-L'-Ut) 

- 

0.25845371 

Selected 

^^rM^.) 

- 

0.153 

t-t-(-fel 

- 

0.10545 

«--«H4:) 

- 

0.22275 

-UU.i^^^) 

- 

0.0357 

-Ap.  (4>x:)p 

rad  -  cm. 

0.2 

Selected 

^./^> 

- 

0.5 

^^,-  C4TD;y 

rad  -  cm. 

0.1 

R^/R 

- 

0.272846 

Ro 

ohms 

1.91  X  10"-^ 

R 

ohms 

7  X  10"-^ 

E 

volts 

min 

0.135  -  max  0. 

225 

I 

Amperes 

min 

15.2   -  max  23 

.6 

w 

Watts 

min 

0.4  4   -  max  1. 

41 

n 

- 

min 

6%    -  max  9% 

TABLE 

10 
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E.   ESTIMATION  OF  AVAILABLE  HEAT  POWER  FROM  THE  EXHAUST  GASES 
The  recently  introduced  catalytic  converter  in  the 

exhaust  system  of  automobiles,  offers  the  advantages  of 

additional  release  of  heat  and  higher  temperatures,  by 

oxidation  of  the  unburned  CO  and  H2 . 

Heat,  also,  can  be  extracted  by  cooling  down  the  gasses 

which  come  out  of  the  output  flange  of  the  converter. 

The  thermodynamic  analysis  which  follows  is  based  on 

3 
data  provided  by  FORD  motor  company.   A  2  300  cm  engine  is 

examined  and  the  results  of  the  exhaust  gas   analysis  before 

and  after  catalytic  converter  are  summarized,  for  three 

different  power  levels,  in  Table  11. 

The  analysis  is  divided  into  three  parts  with  the  neces- 
sary logical  sequence. 

The  first  part  contains  preliminary  calculations,  as 
the  moles/min  and  pressure  of  exhaust  gasses.   The  second  is 
involved  with  the  calculation  of  the  available  heat  energy 
from  the  catalytic  action  of  the  converter,  and  the  third 
is  referred  to  the  estimation  of  the  possible  extracted 
heat  energy  by  cooling  down  the  gasses  after  the  catalytic 
converter. 

1.   Preliminary  Calculations 

a)   From  Table  11  at  idle  speed,  the  engine  produces 
exhaust  gasses  at  a  rate  of  10  ft  /min  or  0.28  32  m  /min. 
One  mole  of  any  gas  at  the  state  of  1  atm  and  2  73  K"  obtains 
a  volume  of  0.0224  m^  and  at  standard  state  -  1  atm,  298  K**  - 
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0.02445  m  .   Therefore,  the  flow  of  the  exhaust  gasses  in 
terms  of  gram-moles  per  minute  is  11.583, 

Based  on  the  percent  mole  analysis  of  Table  11, 
the  actual  amount  of  the  several  constituents  of  the  exhaust 
gas  can  be  summarized  as  follows: 


CURB  IDLE 

60  MPH 

ROAD  LOAD 

PEAK  POWER 

GAS 

gram-mole/min . 

gram-mole/min . 

gram-mole/min 

CO2 

0.80269 

8.2700 

17.316 

CO 

0.4981 

0.3475 

12.787 

H20 

1.0500 

8.548 

24.509 

"2 

0.1865 

0.0 

5.238 

O2 

0.64864 

3.3358 

7.992 

N2 

8.39707 

48.995 

198.472 

TOTAL 

11.583 

69.497 
TABLE  12 

266.404 

I 


b)   The  partial  pressures  of  the  gasses  can  be 
calculated  from: 

^^     "^   Vtot  (5-4  8) 


where : 


M„  =  gram-mole  of  the  gas 
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gas  constant 


T  =  gas  temperature 


TOT 


volume  occupied  by  all  gasses,  comprising 
the  mixture  at  temperature  T. 


With  M^  taken  from  Table  12  the  partial  pressures 
were  calculated  for  three  different  power  levels,  and  the 
results  are  shown  in  tabular  form  as  follows: 


CURB  IDLE 

60  MPH  ROAD  LOAD 

PEAK  POWER 

GAS 

atm 

atm 

atm 

CO2 

0.06931 

0.1190122 

0.066007 

CO 

0.04301 

0.005001 

0.048 

H20 

0.09067 

0.123013 

0.09201 

«2 

0.0161 

0.0 

0.1967 

^2 

0.056007 

0.048005 

0.030003 

N2 

0.725045 

0.70508 

0.7451 

TOTAL 

1.000142 

1.0001112 
TABLE  13 

1.00079 

Based  on  the  data  of  Table  13,  it  is  reasonable  to 
consider  that  the  chemical  process  in  the  catalytic  converter 
takes  place  under  constant  pressure  conditions. 
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2 .   Released  Heat  From  Converter  Action 

The  change  in  enthalpy  produced  by  the  procedure 
of  a  chemical  reaction  from  reactants  in  their  natural  state 
at  25°  C  and  1  atm  to  products  in  their  natural  state  at 
25°  C  and  1  atm  is  referred  to  as  the  "standard  heat  of 
reaction"  at  25°  C. 

The  "standard  heat  of  formation"  is  the  change  in 
enthalpy  resulting  from  the  formation  of  a  compound  at 
25°  C  and  1  atm  from  the  elements  in  their  natural  states 
at  25°  C  and  1  atm.   The  "standard  heat  of  combustion"  is 
the  change  in  enthalpy  resulting  from  the  reaction  of  a 
compound  with  elemental  oxygen  each  initially  at  25°  C  and 
1  atm  to  produce  specifically  defined  products  at  25°  C 
and  1  atm. 

The  heat  of  a  given  reaction  may  be  obtained  by 
constructing  a  set  of  formation  and  combustion  reactions 
such  that  the  sum  of  the  individual  reactions  is  the  reaction 
desired.   The  siom  of  the  enthalpy  changes  of  the  individual 
reactions  will  then  be  the  standard  enthalpy  change  of  the 
given  reaction.   If  the  heats  of  formation  data  are  available 
for  all  compounds  involved,  the  standard  heat  of  reaction 
may  be  obtained  by  subtracting  the  sum  of  the  heats  of  for- 
mation of  the  reactants  from  the  sum  of  the  heats  of  forma- 
tion of  the  products.   Normally,  heats  of  formation  are 
available  for  inorganic  compounds  and  heats  of  combustion 
are  available  for  organic  compounds. 
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Since  the  initial  and  final  states  differ  from  25°  C 
and  1  atm,  the  change  in  enthalpy  may  be  evaluated  by  a 
series  of  steps  so  chosen  that  the  enthalpy  change  of  each 
individual  step  can  determined.   These  steps  are  shown  on 
the  following  diagram. 

The  change  in  enthalpy  for  steps  1  through  3  and 
5  through  16  are  calculated  from: 


^^i'-    ^c/cp-dT 


(5-49) 


where  M.  is  gram-mole  and  c    is  the  molal  heat  capacity  of 
the  gas.   The  molal  heat  capacity  is  a  strong  function  of 
temperature  and  can  be  expressed  as : 


.  =  C<-4ioi:T4C-T^ 


(5-50) 


where  T  is  in  K** . 

The  coefficients  o(.  /  b^  and  c^  are  displayed  in 


Table  14. 

COMPOUND 

FORMULA 

0(i 

b^  X  10 

Carbon  dioxi 

de 

CO2 

6.214 

10.396 

Carbon  monoxide 

CO 

6.420 

1.665 

Water 

H^O 

7.256 

2.298 

Oxygen 

^2 

6.0954 

3.2533 

Hydrogen 

«2 

6.9469 

-  0.1999 

Nitrogen 

N2 

TABLE 

6.4492 
14 

1.4125 

c 

j^  X  10 

- 

3 

545 

- 

0 

196 

0 

283 

- 

1 

.0171 

0 

4808 

- 

0 

0807 

-6 
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step  4  is  calculated  from  tables  where  the  heat  of 
formation  for  different  compounds  is  given.   The  change  in 
enthalpy  is  obtained  by  subtracting  the  heat  of  formation 
of  the  reactants  from  that  of  the  products . 

Table  15  displays  the  change  of  enthalpy  for  each 
step,  for  three  power  conditions: 


CURB  IDLE 

60 

MPH  ROAD  LOAD 

PEAK  POWER 

ENTHALPY 
CHANGE 

Calor/min 

Calor/min 

Calor/min 

«1 

- 

1340.11 

-  1418.0 

-  48025.0 

«2 

- 

955.36 

-   738.7 

-  36487.0 

«3 

- 

492.60 

0.0 

-  20702.0 

«4 

-46303.00 

-21285.0 

-824828.0 

«5 

+ 

2394.40 

+  2114.0 

+  70449.0 

«6 

+ 

728.60 

0.0 

+  25102.0 

"7-12 

+ 

741.00 

-   930.0 

-   6788.0 

«8-13 

+ 

741.30 

-  1157.0 

-   6354.0 

"9-14 

+ 

4859.00 

-  4129.0 

-  44372.0 

"10-15 

+ 

714.00 

-   282.0 

0.0 

"1I-I6 

+ 

3.40 

0.0 

+    885.0 

TOTAL 

-38909.37 

-27826.0 

-892886.0 

TABLE 

15 
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3.   Estimation  of  Available  Heat  from  Hot  Gasses 


I 


The  gasses  flowing  out  of  the  catalytic  converter 
have  a  high  thermal  content,  which  partially  can  be  extracted, 
introducing  a  heat  exhanger  in  series  with  the  catalytic 
converter . 

A  possible  arrangement  of  converter  with  heat 
exhanger  in  series  is  shown  in  Figure  25,   The  change  in 
enthalpy,  because  of  the  cooling  process  of  the  exhaust 
gasses,  can  be  estimated  from  equation  (5-49). 

For  a  temperature  drop  200  F** ,  the  change  in  enthalpies, 
for  different  load  conditions  are  shown  in  Table  16. 


CURB  IDLE   60  MPH  ROAD  LOAD 


ENTHALPY 
DROP 


BTU/min 


BTU/min 


PEAK  POWER 


BTU/min 


CO. 


1737 


10119 


CO 

5 

15 

^^2 

102 

2400 

'^2^ 

1355 

7846 

'"2 

0 

0 

'-2 

6894 

35080 

)TAL 

10090 

55460 

TABLE 

16 

31173 

3468 

0 

30700 

1063 

169892 

236296 


94 


^  - 

U) 

< 

I 

X 

r- 

u            ' 
t-           1 

b)           1 
>           1 
^           1 

<I 

V 

o          1 

y 

FIGURE  25 
CATALYTIC  COMERTER  AND  HEAT  EXCHANGER  COMBINATION 
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In  suiranary,  the  total  amount  of  available  heat  power 
from  the  exhaust  gasses  of  the  examined  engine,  under  the 
assumed  conditions  is  determined  as: 

1)  Idle  condition     38909  +  10090   =   50000  cal/min 

=   3486  watts 

2)  60  MPH  ROAD  LOAD    278264  +  55460  =   83286  cal/min 

=   5808  watts 

3)  PEAK  POWER         892886  +  236296  =  1129182  cal/min 

=  78741  watts 


F.   FINAL  CONSIDERATIONS  ON  GENERATOR 

The  number  of  required  thermoelectric  couples  for  a 
typical  design  of  the  generator  is  based  on  the  available 
thermal  power  which  corresponds  to  60  mph  road  load  of  the 
automobile  engine. 

The  design  specifications  of  the  thermoelectric  generator 
are  shown  in  Table  17. 


96 


ITEM 


UNIT 


VALUE 


REMARKS 


Available  thermal  power 

Efficiency 

Electrical  output  power 

Current 

Voltage 

Number  of  couples 


q>insulation 


watts 

5808 

K» 

900 

SELECTED 

K*» 

335 

SELECTED 

- 

8.9% 

watts 

517 

amps 

22.3 

volts 

23.2 

- 

43.2 

cm. 

8.00 

cm. 

10.36 

rad 

0.2 

rad 

0.1 

cm. 

1.00 

rad 

0.1 

TABLE  17 

The  maximum  working  temperature  of  the  selected  thermo- 
electric materials  is  restricted  to  900  K" .   It  is  possible 
to  have  higher  temperatures  in  the  wall  of  the  catalytic 
converter  when  the  engine  is  overloaded.   In  order  to  protect 
the  hot  junction  thermoelectric  materials  from  damage,  a 
protective  system  must  be  provided.   Such  system  could  be  a 
single  by-pass  line  for  the  exhaust  gasses,  automatically 


97 


operating  when  a  certain  temperature  is  achieved,  or  intro- 
duction of  a  heat  resisting  material  between  the  converter's 
wall  and  the  hot  junctions  of  thermocouples,  which  will 
insure  a  sufficient  temperature  drop. 

The  cold  junction  temperature  was  arbitrarily  selected 
to  be  335  K®.   It  can  be  even  lower,  depending  on  the  heat 
rejection  system  which  is  employed,  i.e.  water  cooling,  forced 
or  free  air  convection. 

The  hot  junction  temperature  is  varying  according  to  the 
developed  power  level  of  the  engine.   This  results  in  the 
variation  of  all  design  specifications  of  the  thermoelectric 
generator,  as  delivered  useful  power,  voltage,  current,  etc. 
For  this  reason,  an  electrical  power  conditioning  device 
becomes  necessary.   The  function  of  this  device  depends  on 
the  special  application  of  the  generator.   That  is,  if  the 
generator  is  employed  as  the  electrical  power  source  of  the 
automobile,  constant  voltage  and  variable  current  is  required, 
and  when  the  generator  supplies  part  of  the  electrical  power 
to  the  thermoelectric  cooler,  constant  current  and  variable 
voltage  is  desired. 

Contact  resistance  and  conductor  losses  were  not  taken 
into  account,  because  the  estimation  of  the  former  is  a  matter 
of  manufacturing  procedure  and  the  latter  depends  on  the 
proper  material  for  the  applied  contact  technique. 

Nevertheless,  a  reasonable  figure  of  contact  resistance 

2 

is  of  the  order  of  10  micro-ohms  per  cm  ,  which  is  very  small 

compared  to  the  resistivity  of  the  thermoelectric  materials. 
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Figure  26  shows  the  arrangement  of  catalytic  converter  with 
the  attached  thermoelectric  generator. 
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FIGURE  26 
T.E.  GENERATOR  ASSEMBLY 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  obtained  figures  for  the  thermoelectric  cooler  and 
generator  could  be  considered  as  quite  promising  at  least 
from  the  standpoint  of  feasibility  of  the  application. 

The  thermoelectric  cooler  for  7000  BTU/hour  was  found 
to  require  19  53.0  watts  of  electrical  power.   From  this  figure 
50  watts  are  dissipated  for  the  required  air  blowers.   If  the 
four  units,  which  comprise  the  device,  are  fed  in  parallel, 
then  an  external  conventional  generator,  coupled  to  the 
engine,  is  needed  with  33.  volts  and  57.6  amperes  as  rating 
characteristics . 

When  the  thermoelectric  generator  and  the  cooler  are 
considered  as  one  system,  then  the  required  external  electrical 
power  is  reduced  by  18  to  25%,  depending  on  the  operating 
power  level  of  the  engine,  and  correspondingly  the  coefficient 
of  performance  of  the  cooler  is  increasing. 

A  direct  comparison  between  the  conventional  air  condi- 
tioning system  and  the  thermoelectric  cooler-generator  scheme 
is  not  easy,  because  there  are  not  available  enough  data  for 
the  latter.   In  Table  18  an  attempt  for  comparison  is  made  in 
the  most  broad  characteristics. 
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CHARACTERISTIC 


CONVENTIONAL 
SYSTEM 


THERMOELECTRIC 
SYSTEM 


Voliome 

Weight 

Complexity 

Maintenance 

Noise  level 

COP 

Working  fluids 


Both  systems  have  almost  the  same 
worst  best 

complex  simple 

required  almost  none 

relatively  noisy  quiet 

Both  systems  are  in  the  same  level 
yes  no 


TABLE  18 


I 


From  the  standpoint  of  manufacturing  cost,  the  existing 
figures  show  that  the  thermoelectric  system  is  two  to  three 
times  more  expensive  from  the  corresponding  conventional. 
It  must  be  mentioned  that  the  cost  of  the  conventional  system 
is  referred  to  the  "series  production  cost"  and  the  corre- 
sponding for  the  thermoelectric  system,  to  experimental 
constructions.   It  is  expected  that  whenever  the  thermoelectric 
system  will  be  manufactured  in  series,  the  cost  will  decrease 
enough  and  eventually  it  will  be  compatible  or  even  cheaper 
from  the  conventional  system,  since  the  thermoelectric  material 
itself  is  not  expensive. 

The  analysis  showed  that  the  efficiency  of  the  generator 
depends  strongly  on  the  thermoelectric  material.   The  obtained 
figures  —  6%  minimum  and  9%  maximum  —  for  the  generator's 
efficiency,  are  considered  satisfactory  compared  to 
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efficiencies  reported  already  in  the  literature.   Bates  and 
Weinstein  developed  a  similar  design  using  PbTe  and  SiGe 
which  complemented  each  other  over  the  temperature  interval 
from  50C°    to  1000  C° .   Their  generator  achieved  an  efficiency 
of  8.7  percent  with  a  hot  junction  temperature  of  800  C® 
(1073  K°) .   The  present  design  shows  an  efficiency  of  9% 
with  Tjj  =  940  K°  and  T^  =  335  K°. 

The  requirements  for  the  primary  electrical  source  of  a 
small  automobile  vehicle  are  satisfactorily  fulfilled  from 
the  specifications  of  the  segmented  type  thermoelectric 
generator. 

Table  19  displays  the  electrical  power  requirements  for 
a  small  size  car. 


ITEM 

WATTAGE 
(12-Volts) 

high  beam 
Head  lamps 

low  beam 

45 
40 

Front  lamps 

21 

Tail  lamps 

21 

Other  lamps 

10 

Wiper  motor 

10 

Electrofan  motor 

11.5 

TOTAL  158.5 


TABLE  19 
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The  thermoelectric  generator  provides  2  30  watts  in  the 
idle  condition  of  the  engine  and  517  watts  in  the  60  mph  road 
load.   So,  from  this  simple  examination,  the  thermoelectric 
generator  can  replace  the  conventional  alternator,  offering 
the  great  advantages  of  fuel  economy  and  elimination  of 
moving  parts. 

The  realization  of  the  discussed  applications  requires 
further  theoretical  and  practical  work  to  be  done.   The 
systems  must  be  examined  thoroughly  from  the  following 
standpoints : 

1)  Stress  analysis  of  T.E.  couples  for  T.E.  generator 
and  T.E.  cooler. 

2)  Thermal  and  electrical  insulation. 

3)  Contact  techniques. 

4)  Fabrication  of  thermoelectric  pellets. 

5)  Power  conditioning  circuits  for  T.E.  generator. 

6)  Control  system  for  the  automatic  co-operation  of  cooler, 
T.E.  generator  and  conventional  alternator. 

It  must  be  stressed  that  the  present  work  was  developed 

under  the  scheme  of  the  feasibility  study  instead  of  a 

detailed  design. 
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APPENDIX  A 
APPROXIMATE  ANALYTIC  SOLUTION  OF  HEAT  BALANCE  EQUATION 

The  heat  balance  equation  is: 

^     dr^   dT     \dt~J  ^^  J     ^r  dr     ^r 

where:  ^^C^"^  r   fCT)  /  €(t)  f      are  functions  of  temperature. 
Equation  (A-1)  can  be  written  also  as: 

XT)^V-i!l  .  M.  W^-(dl.j,  ^CT^^VJI  1  l6CT)-mril  ,fp(T>0(A-2) 
dr"^   oT      ^dr '        dr        dr   ^ 


Expression  T -r—       can  be  written  as  —7-      or 


If  the  variable  x  is  introduced  such  that  x  =  In  r   then: 


^  dir  ^  dr. 
dr   "    dx 

dT   _  1  dr 
dr  ~   '^  cAx 


(A-3) 


(A-4) 


drLdr]      dv-LrdvJ     dr»-rleix        r  d^x^    dr 


105 


cJ\  .   1  Jt  .  1  J't 


dlr-^-     r^dx   r^  dx^  **'^' 


Applying  transformation  (A-7)  to  equation  (A-2) ,  it  becomes: 


or; 


(A-8) 


Dividing  by  the  product  ml  equation  (A-8)  becomes: 


The  following  transformation  is  introduced: 
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Differentiation  with  respect  to  T,  equation  (A-10)  becomes 
or 


Applying  the  transformation  (A-11)  to  equation  (A-9) ,  the 
latter  becomes: 


i^.  ecTncTKsCT^=.o  (a-i2) 

ci-T      :z- 


Through  the  attempted  transformation  the  order  of  the  original 
differential  equation  was  reduced  but  the  non-linearity 
remained,  so  the  solution  of  it  at  least  is  not  obvious.   From 
this  point,  the  analysis  of  the  problem  is  based  on  the 
assumption  that  in  a  first  approximation  the  electric  effects 
are  independent.   That  means,  it  is  assumed  that  the  Thomson 
effect  is  not  interacting  with  the  Joule  effect. 

From  this  standpoint,  the  electric  effects  can  be  studied 
separately  and  their  influence  can  be  added  finally  in  the 
interpretation  of  the  formulae  for  the  rate  of  heat  flow  in 
the  two  ends  of  the  semiconductor  element. 
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1.   Influence  of  Thomson  Effect  on  Heat  Conduction 

It  is  assumed  that  the  Joule  effect  is  not  present  in 
the  system.   Then  it  can  be  set  P(T)  =0. 
Equation  (A-12)  becomes,  for  one  possible  current  direction: 


^^  -^  6CT)=  0  (A-13) 

elT 


Integrating  and  introducing  the  transformation  (A-10) , 
equation  (A-13)  becomes; 

-^CT)nnn^:.B-l/^CT)clT  (A-14) 

where  B  is  the  integration  constant  which  must  be  determined. 
Transforming  the  independent  variable  x  to  the  original  r, 
equation  (A-14)  becomes: 

■>^(J)Tt^^iI  ^  b-l/6CT)cJT  (A-15) 

The  boundary  conditions  of  the  system  are: 


r  =  R^  T  =  Tjj 


r  =  R2  "^  =  ^c 


Under  these  conditions,  the  rate  of  heat  flow  in  the  two  ends 
of  the  element  is  given  by: 
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(A-16) 


-^^-^^'-^^l.u.-'^ 


Equation  (A-15)  can  also  be  written  as: 


^eT:)'mril  ^  B  (i-  i  /Vc-DdT  )  (a-i7) 


Integrating  both  sides: 


(A-18) 


The  following  notations  are  used: 


Q  =   rate  of  conduction  heat  flow  = 


lAf) 


i   '^McJt 


v.. 


q   =   rate  of  Thomson  heat  flow 


From  equations  (A-16)  it  can  be  seen  that  B;^Q  . 

The  equality  sign  holds  for  the  low  temperature  end. 
Also,  for  normal  conditions,   Q>  |  q„  \     and  consequently, 
— 5 —  ^  1  ,   which  justifies  the  required  conditions  in  order 
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to  expand  the  fraction   — in  a  series  of  the 

following  form: 


i-riT-:;:;  ~-  '^Tt'^'V  lii^-^^f^-^if^^-^^ir-is: 


Equation  (A-18)  now  can  be  written  as: 
^  =  fe  i  lCT)[l.(-£-/lcT)dT)4  (yfCr)dTjV..^(-j/c6)dT)  L^(A-20) 


Integrating  term  by  term  and  after  some  algebra,  equation 
(A-20)  can  be  written  as: 


where: 


for  i  =  1,2, ... ,n 

From  equation  (A-20)  it  can  be  set,  for  a  first  approximation, 
B  =  Q.   Then,  B  for  the  second  approximation  becomes: 
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(A-22) 


It  is  sufficient  to  keep  only  the  first  two  terms  of  equation 
(A-22),  since  higher  order  terms  are  quite  small. 


B 


^n.-^ 


6CT)c(Td- 


y   ':\ct)cIt 


(A-23) 


Substitution  of  (A-23)  in  equation  (A-16)  yields: 
'XT^I-mr^    ^—^        ')CT)dT  -  i  ^^1^= lis 


(A-24) 


/:)"ct)J- 


2.   Influence  of  the  Joule  Effect  on  Heat  Conduction 
Here  it  is  assumed  that  the  Thomson  effect  is  not 
present,  so  setting  &(T)  =0   equation  (A-12)  becomes: 


^    ecT)UT) 


=  o 


(A-25) 
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Rearranging  and  integrating,  it  is  obtained: 

T 

Z'-^  D''-2y7^CT)(>CT")ciT 

where  D  is   the  integration  constant. 

Introduction  of   the   transformations: 


2    =   •— 


dx 


in  the  above  equation  yields: 

The  appearance  of  two  signs  in  equation  (A-2  6)  leads  to  the 
conclusion  that  the  temperature  distribution  along  the  element 
will  present  a  local  maximum  T  at  a  point  corresponding  to 
some  radial  distance  R  . 

Rearranging  equation  (A-26) ,  integrating  both  sides 
and  applying  the  boundary  conditions,  the  following  expression 
for  D^is  obtained: 


D,  -  ;^^  1 1  r    }^]t^ ;„  -  f^T^. -..,  I  '^-2^> 


^'Ml)  Wv-%^h^M'^  ([-|i;eco;,c.)df . 
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and  for  the  point   (R_,  T  ) 
mm 


D,=  l[2jr'^tT^^CTWTT'"  (a.28) 


There  is  a  certain  value  of  the  current  which  could 
establish  the  condition  of  no  back  heat  flow  in  the  hot  end 
of  the  element.   That  means,  it  is  possible  to  have  the  point 
(Rj^,  T^)  shifted  to  coincide  with  the  point  (P,  ,  T  )  ,  applying 
a  certain  value  of  current. 

This  value  of  current  can  be  found  from  equations 
(A-27)  and  (A-2  8)  and  it  is: 


-^     /    :\Ct)cIt 

(A-29) 


^<')  k[^  C^C.,,c.,A-rr 


The  above  result  might  be  important  for  the  thermo- 
electric cooler,  which  should  operate  under  the  most  efficient 
conditions. 

The  rate  of  heat  flow  at  the  two  element  ends  is  given 
by  the  following  equations. 


(A-30) 


r^\        =  Q 


1 

2. 


113 


It  is  obvious  from  equation  (A-26) ,  in  order  for  the 
left  side  to  be  real,  the  following  condition  should  exist: 

D^>  OJ?  f   (^CT^^^CT-)  dT  (A-31) 

for  the  range  of  temperature  under  consideration. 
Also,  from  equation  (A-26)  by  integration  of  both  sides 


r>  -   ^    \      I  ^Ct^cIt     1  (A-32) 


'  ^^(^)Kb-f/;ec-:.^w^.Y 


The  condition  declared  by  equation  (A-31)  justifies 
the  series  expansion  of  the  integrand  of  equation  (A-32) ,  in 
the  binomial  form. 

{^i-.2:^  f  ^iT^'\L'^^^-A        -  C^l-^V)''  (A-33) 

where  p  =>|  and  x  =  -  2^  /  (^ct^) ->^CT)  ciT 

2.  I  ^  I 
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Application  of  the  above  expansion  in  equation  (A-32)  and 
integration  of  it  term  by  term  would  yield,  after  algebraic 
manipulation,  the  following  expression  for  D: 


where : 


(A-34) 


om 


-      0  f^rJ^  ^^^^^^ 

C    =   ^-»vCa^-0  (-m-a-)  ■■■  (l-rx-(.-n-,)\ 


'Jn 


111, 


Following  the  same  procedure  as  before,  one  can  set 
D^=  Q  as  a  first  approximation.   Applying  this  value  of  D  in 
equation  (A-34),  D^ can  be  calculated  in  a  second  approximation 


t)^  =  Q  [  1  -V  -lik_  _, 


3  X. 


%.      ^^ 


(A-35) 


The  higher  order  correction  in  equation  {A-35)  might  be 
considered  to  be  quite  small  to  be  taken  into  account  in  an 
approximation  solution. 
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It  is  sufficient  to  approximate  D  by: 

i""  ^  ""  %ji  (A-36) 

Applying  the  above  expression  for  D^  in  equations  (A-30)  the 
rate  of  heat  flow  at  the  two  ends  of  the  element,  now  it  is 
given  by: 


(A-37) 


I 


3.   Composite  Solution 

Based  on  the  assumption  that  the  electric  effects  — 
Joule  and  Thomson  —  are  non-interacting,  superposition  of 
the  obtained  results  yield  the  following  composite  solution 
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for  the  rate  of  flow  of  conduction  heat  in  the  semiconducting 
pellet: 


-^(T) 


TTt' 


dr 


(A-38) 


^L-n^r^ 


TM   ^T 


-m 


^C 


T.   Act)  /.eCT)ofToI-^ 
_  i>^is! lis 


(A-39) 


Equation  (A-38)  is  quite  complicated  to  be  used  in  the  analysi; 
of  the  thermoelectric  cooler  or  generator,  and  it  can  be 
replaced  by  a  simpler  but  less  exact  expression  as: 


I  6CT)/:|CT)dTclT 


where : 


hi 


(A-40) 
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The  introduced  error  in  equation  (A-40)  has  been 
investigated  using  computer  techniques  and  was  found  to  be  in 
the  order  of  +5.8%  and  -3%  for  the  resistivity  terms  of  (A-40) 
for  P  and  N  type  semiconducting  pellets. 

Analytically,  a  semiconducting  SiGe  allow  has  been 
used  with  electrical  resistivity  and  thermal  conductivity  as 
shown  in  Figures  2  7  and  28.   For  T„  =  645  K"   and  T^^  =  400  K° 
the  following  results  were  obtained: 

P-TYPE  N-TYPE 


P  =  1.12535  X  10  -^      f  =  0.64771  x  lO""^ 
^  =  1.275  X  10"^      if"  °*^^  ^  ^°~^ 


where: 


These  small  errors  have  a  very  slight  effect  in  the 
calculated  amount  of  the  conduction  heat  rate  of  flow  as  it 
is  calculated  by  (A-38)  and  (A-40) .   The  calculated  error 
was  of  the  order  of  0.2%. 

The  validity  of  the  assumption  of  non-interacting 
electrical  phenomena,  in  a  first  approximation,  can  be 
justified  by  replacing  the  arbitrary  relation  pCT  )  ^ '^CT")  o<  CT) 
with  constant  values  in  equations  (A-38)  and  (A-39) .   Such  a 
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FIGURE  27 

THERI.IAL  CONDUCTIVITY  x  10^  -  WATT/CM/DEG. 


119 


.     . 

/ 

t^-i^p^y 

^ 

/ 

X 

y 

/. 

/ 

^^^ 

p-TYP^,^-'''^ 

'^ 

200 


400       600       800      1000      1200 


Temperature  K 


FIGURE  28 


ELECTRICAL  RESISTIVITY  -  OHM  -  CM  x  10- 
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substitution  gives  the  exact  form  of  the  rigorous  solution 
of  the  simplified  differential  equation  as  obtained  in 
Appendix  B. 
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APPENDIX  B 
SOLUTION  OF  SIMPLIFIED  HEAT  BALANCE  EQUATION 

The  heat  balance  equation  in  the  thermoelectric  element 
with  constant  properties  was  defined  in  Section  II  as: 


£l  ^_Lil  ^_B_  ^ 


O  (B-1) 


dr'^   r  dr    V" 

n     ei!_ 

where  B  =  -^- -_ 


2 

Multiplying  by  r  ,  equation  (B-1)  becomes 


S-'^Al  4- r-^  4-B=  O  (B-2) 

Introducing   the   transformation      x  =   In  r    ,      equation    (B-2) 
becomes : 

£l    -  _B  (B-3) 

Integration  of  equation  (B-3)  yields: 

^I  ^  _'^x  -+  C^            ..-  (B-4) 

T=  C^-+C,X  -Nl^X"    _-  -                 (B-5) 
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Replacing  x  with  In  r  equation  (B-5)  becomes: 


T=  C^  -V  q  (r^r  ~  -^^  {^'^^^  (B-6) 


The  initial  conditions  for  the  system  are; 


r  =  R^  T  =  Tg 


r  =  R2  T  =  Tj, 


Introducing  the  above  conditions  to  equation  (B-6) ,  it  yields: 
Th  =C:,^C,I^R,  -lB(l.rvC-,)^      -        (B-7) 
■Tc-C^-fcJ-r,U;,-lB(^'flRiT  (B-8) 

From  equations  (B-7)  and  (B-8)  the  constants  of  integration  C-, 
and  C2  are  determined  as : 


C.  - 


-.       zL  -•■-  -  -'    ^-^^    ^WJ   -•  __(B_9) 


c^ .  Th  - c,  -L^^  -V  :|  [Cco]^  (B-io) 
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The  rates  of  heat  flow  at  the  hot  junction  (T  =  T„)  and  at 

n 

the  cold  junction  are  defined  as: 


(B-ll) 


(B-12) 


The  temperature  distribution  along  a  SiGe  segment  is  shown 
in  Figure  29. 
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